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1. INTRODUCTION 

The study of jet production in association with 
a vector boson at hadron colliders has been suc- 
cesful in the recent past. The advantage of this 
signal over pure jet production is that the lep- 
ton(s) from the vector boson decay can be used 
as a trigger, such that jets can be studied free 
from jet-trigger bias. Furthermore, the lower rate 
obviates the need for prescaling. On the theory 
side, progress in calculational techniques to con- 
struct next-to-leading order (NLO) Monte-Carlo 
programs Q has allowed a meaningful confronta- 
tion with data pi . For a review of such techniques 
and more recent progress, see ||. 

The tagging of heavy hadrons in the jet offers a 
unique possibility of studying the hadronic struc- 
ture inside the jet. By considering jets where 
the leading hadron is tagged, a clear connection 
can be made with perturbative QCD. At the par- 
ton level, the tagging of a heavy hadron corre- 
sponds to the tagging of a heavy-flavor quark. 
Experimentally, the presence of a D or B meson 
is inferred through its decay products. The im- 
portance of heavy flavor tagging has been clearly 
demonstrated in the analysis that led to the top 
quark discovery jl. In the future, heavy flavor 
tagging will continue to be an important analysis 



tool. It will provide more detailed information 
about the event, and test of the underlying QCD 
theory. 

If one demands the presence of a charm quark 
in the jet recoiling against a W-bosou, the signal 
is directly sensitive to the strange quark distribu- 
tion function in the proton, at a scale of the order 
of the W-mass. 

Replacing the charm quark by a bottom quark, 
one could use this process as an alternative cali- 
bration of the 6-quark tagging efficiency. This will 
be useful at luminosities achieved by the Main 
Injector. However, by far the dominant contribu- 
tion to VF+bottom production is due to W + bb 
production, where the heavy quark pair is pro- 
duced by gluon splitting. Here the inclusion of the 
gluon to B meson fragmentation function is prob- 
ably more important than the inclusion of NLO 
effects to the ly+bottom process, which implies 
that this reaction can constrain this fragmenta- 
tion function. 

In this report, which is based on Ref. JtJ, we 
present the calculation of the QCD corrections 
up to 0{a 2 s ) of the process pp — > W + Q where Q 
is an heavy quark, keeping the mass of the heavy 
quark explicit, and the final state fully exclusive. 
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2. METHOD 

The method we used to calculate the 0(a 2 s ) QCD 
corrections to pp — ► W + Q is a generalization Q 
of the phase space slicing method of Ref. to 
include massive quarks. 

The leading order (LO) calculation involves the 
simple subprocess q + g — > VF + Q, where g is a 
light quark and g a gluon, at tree level. 

The virtual corrections consist of the inter- 
ference between the lowest order diagrams and 
their one-loop corrections. The integrals over the 
loop momenta were performed in d = 4 — 2e di- 
mensions. Ultraviolet singularities were absorbed 
through mass (on-shell scheme) and coupling- 
constant (modified MS scheme |10|]) renormaliza- 
tion. The remaining soft and collinear singular- 
ities, appearing as 1/e 2 and 1/e poles, factorize 
into a universal factor multiplying the Born cross 
section. 

The real corrections consist of the contributions 
from all the subprocesses ij — > W Q k, where 
k are massless partons, and the subprocess 
ij — » W QQ plj . Some of these contributions 
exhibit soft and/or collinear singularities. We do 
not consider diagrams where the heavy quark is 
in the initial state. The method we used to isolate 
the singularities consists of slicing up the phase 
space, using a cut-off s m in and color-ordered sub- 
amplitudes, into a hard region (s m i n < 2p r -p s for 
any r, s), containing no singularities, and a com- 
plementary region in which the final state parton 
is either soft or emitted collinearly with one of 
the initial state partons. In the hard phase space 
region, one can work in four dimensions and per- 
form the phase space integration numerically. In 
the soft and collinear region, the integration is 
done analytically in d dimensions using soft and 
collinear approximations, which are valid in the 
limit that s m j n is small. The cross section in this 
region again factorizes into a universal factor mul- 
tiplying the Born cross section. The initial state 
collinear singularities are factorized into parton 
distribution functions in the MS scheme, using 
the formalism of crossing functions |Q. 

Note that the process ij — ► WQQ is quite 
different from the other subprocesses: the heavy 
quark does not originate from the W vertex, and 



it is independent of s m ; n , because it is free from 
singularities. 

Adding the real and virtual corrections leads to 
the cancellation of all remaining singularities. We 
checked gauge invariance for both the virtual and 
real corrections. Finally, we constructed a Monte 
Carlo program for the present process, including 
these corrections. 

Before showing any numerical results, we first 
list here the default choices we made for param- 
eters and cuts in producing the results of this 
paper. Any deviation from these choices will be 
indicated explicitly. For the case of charm (bot- 
tom) we assumed three (four) light flavors. We 
used both at LO and NLO the CTEQ3M @ 
set of parton distribution functions, and a two- 
loop running coupling constant with Aqq 4 ^ the 
value supplied with the CTEQ3M set (0.239 
GeV). We implemented continuity across heavy 
flavor thresholds |l3| using the parametrization of 
Ref. jnj. We used the Snowmass convention ]l5| ] 
for the definition of a jet. Our conditions on the 
transverse energy and pseudorapidity of the jet 
were £r(jet) > 10 GeV, \rij et \ < 3, and no cuts 
on the W. We took a jet cone size of AR = 0.7, 
the W mass equal to mw = 80.23 GeV, and the 
heavy quark mass m equal to 1.7 GeV for charm 
and 5 GeV for bottom. We used V cs = 0.97 and 
V c d = 0.22 for the relevant Cabibbo-Kobayashi- 
Maskawa matrix elements. We chose the factor- 
ization scale equal to the renormalization scale 
and denote it by fj,, taking [i = my/. At least one 
heavy quark was required to be inside of the jet, 
with the sign of its electric charge correlated with 
the W charge, as in the LO case. 

We verified the independence of the calcula- 
tion on the arbitrary parameter s m i n . The results 
given in the remainder of this report are averaged 
over s m i n between 1 and 10 GeV. 

We further found the scale dependence of the 
inclusive cross section to be somewhat reduced by 
the inclusion of the NLO corrections. 

3. STRANGE QUARK DISTRIBUTION 
FUNCTION 

We first discuss the effect that the inclusion of the 
NLO QCD corrections to W + charm-tagged jet 
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Table 1 

The W + charm-tagged one-jet inclusive cross section in pb for LO, W + QQ, and NLO (including the 
W + QQ contribution) using different sets of parton distribution functions. The statistical uncertainty 



from the Monte-Carlo integration is less than 1%. 

set mass (GeV) LO WQQ NLO 

CTEQ1M m c =1.7 96 20 161 

MRSDO' m c =1.7 81 20 138 

CTEQ3M m c =1.7 83 20 141 

CTEQ3M m b =5.0 0.17 9.09 9.33 



production has on constraining the strange quark 
distribution function s(x,n) in the proton. Here 
x is the momentum fraction of the strange quark 
in the proton, and fi is the factorization scale. See 
Ref. || for an extensive study of this issue using 
the shower Monte-Carlo program PYTHIA || . 

In Table 1 we give the NLO cross section for 
the parton distribution function sets CTEQ1M 
16) and MRSDO' @. The MRSDO' set derives 
its strange quark distribution from di-muon data, 
whereas the CTEQ1M set uses DIS data. At low 
/j, 2 and low x the difference is as much as a fac- 
tor of two. Also shown is the result obtained 
with the more recent CTEQ3M set, which uses 
the same assumption about the strange quark dis- 
tribution as MRSDO'. Comparing the CTEQ3M 
and MRSDO' sets, we see that the difference due 
to using more recent data sets in the global fit for 
the parton distribution functions is small. This 
is also reflected in the cross section for W + cc, 
which is the same for all three sets. We can con- 
clude that the difference between CTEQ1M and 
MRSDO' of 15.4 % is due to the strange quark 
distribution function. This difference becomes 
14.5% when one includes the W + b background 
(9 pb, almost all of it coming from the gluon split- 
ting contribution, see Table 1), and assuming con- 
servatively that each bottom quark is mistagged 
as a charm quark. This shows that the conclu- 
sions reached in Ref. || are still valid at NLO. 
In both the NLO calculation and the PYTHIA 
analysis about 50% of the contributions are ini- 
tiated by strange quarks. One major difference 
is that PYTHIA suggests that the gluon splitting 
contributes about 35%, whereas in the NLO cal- 
culation it is only about 15%. This number is 
however quite sensitive to the choice of scale jj, in 



the gluon splitting contribution. 

4. MASS DEPENDENCE OF JET PRO- 
DUCTION 

Next we compare W + untagged jet production^ 
and W + charm-tagged jet production. In 
particular, we study the jet transverse energy 
(£r(jet)) distribution. In Fig. |l| we show the ra- 
tio of the charm-tagged jet over the untagged jet 
-Er(jet)-distribution for the LO and NLO cases. 
At LO the charm-tagged jet is simply represented 
by a charm quark. The ratio in Fig. [l] has a char- 
acteristic shape which can readily be understood 
at tree level. At low Et (jet) the charm-tagged jet 
rate is suppressed relative to the untagged-jet rate 
due to its fermionic final state. The untagged- 
jet rate is dominated by the gluonic final state 
which has a soft singularity, that is absent for 
the fermionic final state. At high Et (jet) we 
again observe a relative suppression of the charm- 
tagged jet because at LO this process has a gluon 
in the initial state. At high Et (jet) the dom- 
inant scattering in the untagged-jet rate is due 
to quark- antiquark collisions, again favoring the 
gluonic final state. Clearly, apart from an ap- 
proximate overall if-factor, the NLO cross sec- 
tion retains these features. Also shown in Fig. [j] 
is the W + cc contribution. At low Et (jet) its 
suppression is more pronounced due to the charm 
quark pair production threshold. At high Et (jet) 
there is no suppression for this process because it 
is instigated by a quark- antiquark collision. 

We now turn to the mass dependence of the 
cross section. We will show that there are im- 



2 For the untagged process we take five massless quark 
flavors. 
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Figure 1. Ratio of W + charm-tagged inclusive 
one-jet production over W + inclusive untagged 
one-jet production as function of the jet trans- 
verse energy. The solid line is the NLO ratio and 
the dashed line the LO ratio. The W + cc con- 
tribution to the NLO ratio is also shown (dotted 
line). 



portant mass effects at NLO, especially when we 
look in detail at the tagged jets. (At LO one may 
safely take the mass to zero; the answer does not 
change.) 

Let us first consider the W + QQ contribution. 
Because the mass of the heavy quark regulates 
the collinear singularity, it is expected that the 
strongest mass dependence will come from the 
collinear region. In this region the cross section 
factorizes into the cross section for gluon pro- 
duction multiplied by a universal factor. After 
integration over the invariant mass of the heavy 
quark pair we find that the mass dependent part 
of this universal factor has the following form: 



(z)ln( ^4 I 'I: 



(1) 



where M is the upper limit of the heavy quark 
pair invariant mass defining the collinear region, 
and P q q^g(z) is the massless Altarelli-Parisi Q 
splitting function: 



We chose the following definition of z 
E + P n 



Ejet + Pje 



(3) 



where E and Pm are the heavy quark energy and 
momentum projected on the jet direction, and 
Ej et and Pj e t are the jet energy and momentum. 
Other choices, such as z — By(Q)/£^r(jet), do 
not change any of the conclusions in what fol- 
lows. In the strictly collinear limit M is much 
smaller than the energy of the gluon, but in a 
leading logarithmic approximation one may take 




(2) 



Figure 2. a) Ratio of the W + cc component of 
the charm-tagged one-jet inclusive cross section 
to the W + gluon cross section, as a function of 
the jet tranverse energy, b) The z-distribution of 
the W + cc component. 



M to be of the order of £/r(jet). The behavior 
of Eq. (|l|) can be seen qualitatively in Fig. ||a, 
where the ratio of the W + cc cross section over 
the W + gluon cross section is shown as a function 
of the transverse energy of the jet. One can see an 
approximate logarithmic enhancement with in- 
creasing Et, as predicted by the leading logarith- 
mic approximation. However the z-distribution, 
plotted in Fig. ||b, does not conform with the in- 
dependence in Eq. ([!]). First, the peak at z = 1 is 
due to events where the Q is not inside of the jet, 
such that the whole jet is formed by the lone Q. 
Second, the cross section is suppressed near z = 
and z — 1 (excluding the peak) due to terms in 
the collinear region that depend strongly on z, 
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but not on M. One way to enhance the effect of 
the leading logarithm term of Eq. ([!]) in the z dis- 
tribution is to lower the mass of the heavy quark 
in our calculation. This is done in Fig. ||a, where 
we show the z-distribution for the case m — 0.01 
GeV. Clearly it now resembles the functional form 
of Eq. (|J) much closer. 

The ln(m 2 ) term in Eq. (|l|) diverges in the limit 
of vanishing quark mass, and is not cancelled by 
any other contribution ^. In principle, any ob- 
servable should be "collinear safe", i.e. if the mass 
is taken to zero, the observable should be finite 
and approach the massless result. This is needed 
to describe situations where the relevant scale is 
much larger than the heavy quark mass. In the 
present case we are dealing with a final state di- 
vergence in ln(m 2 ) that should be factorized into 
the fragmentation function of a gluon into a heavy 
hadron. The evolution of fragmentation functions 
will then resum the large logarithms In(E^). This 
problem is the final state version of the problem 
of heavy quark distribution functions |Te[ | . It will 
be discussed in more detail in [§J. Some studies 
in this regard were done in Ref. [p5| . Here, we 
simply keep the mass finite. Let us now turn to 
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Figure 3. The z-distribution of the charm-tagged 
one-jet inclusive cross section, with m = 0.01 
GeV. a) Wcc component, b) Total contribution 
minus the Wcc component. 



3 In the W +1 jet calculation, this singularity is cancelled 
by a companion collinear singularity in the quark-loop cor- 
rection to the outgoing gluon in the Born diagram. In our 
calculation this diagram is not present. 



the NLO single charm contribution (excluding the 
W + cc contribution). In the collinear region, the 
cross section again factorizes, leading, after inte- 
gration over the invariant mass of the collinear 
partons, to the universal factor: 



M 2 



dz 



(4) 



where P qg ^ q (z) is the splitting function: 



Pa. 



Urn 2(1 



l + z 1 



9(1- z- 6) + (| + 2bxS) 5(1 - zj) . 



(5) 



In Fig. 6b we show for the single charm quark 
contribution the z-distribution, with the contri- 
bution of Eq. (^) enhanced by taking m = 0.01 
GeV. Note that it resembles the functional form 
m Eq. (§). From Eq. (||) one derives 



P qg ^ q (z)dz = , 



(6) 



which must hold for the probability to find a 
quark in a quark of the same flavor to be one Jl8[ . 
From Eqs. (Q) and (^) we can now make the im- 
portant observation that as long as the cuts on the 
heavy quark-tagged jet are such that all z-values 
are allowed to contribute, there are no large loga- 
rithms \n(E^/m 2 ). However, if the cuts are such 
that the z-integration is restricted, or convoluted 
with a z-dependent function, some lu.(E^,/m 2 ) 
terms will remain. An example of such a convo- 
lution is the Et distribution of the heavy quark 
itself. All this is illustrated in Fig. |]a, where 
two cases of z-restrictions (z > 0.9 and z < 0.9) 
are plotted in addition to the all-z case. Note 
that for the all-z case the ratio is indeed essen- 
tially independent of Er(jet), but that for the 
z-restricted cases the logarithmic dependence is 
apparent. For completeness we show in Fig. [|b 
the z-distribution for the single charm contribu- 
tion for m = 1.7 GeV. 

Thus care must be taken when calculating 
tagged cross sections in determining whether or 
not there are large logarithms present due to re- 
strictions on the z-integration. Such restrictions 
would also necessitate the introduction of the ap- 
propriate fragmentation function to absorb the 
ln(m 2 ) terms, as in the WQQ case. 
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Figure 4. a) Et distribution ratio of NLO 
charm-tagged one-jet inclusive cross section to 
LO one. Solid line has no restrictions on the z- 
intcgration, the dotted and dashed line have re- 
strictions of z < 0.9 and z > 0.9 respectively, b) 
z-distribution of NLO charm-tagged one-jet in- 
clusive cross section. The W + cc component is 
not included in these plots. 



Note that for the example of the jet transverse 
energy distribution given in Fig. 4 there is no con- 
straint on the z-integration, so that the only final 
state ln(m 2 ) term arises from the Wcc contribu- 
tion. 

5. CONCLUSIONS 

We have presented the calculation of the 0(a s ) 
corrections to the reaction pp — + W + Q. We 
found that the inclusion of the NLO corrections 
does not change the conclusions of Ref . Q about 
constraining the strange quark distribution func- 
tion using W + charm-tagged jet events at the 
Tevatron. However, using our NLO calculation, 
this procedure can constrain the NLO strange 
quark distribution function, once a reasonable 
data sample is collected. We studied the Et dis- 
tribution of the jet containing the heavy quark 
and the mass dependence of the cross section, and 
noted the need to include the heavy hadron frag- 
mentation functions. 
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